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USER’S MANUAL FOR FOMOCO UTILITIES — 

FORCE AND MOMENT COMPUTATION TOOLS FOR OVERSET GRIDS 

William M. Chan * and Pieter G. Buning 
Ames Research Center 


Summary 

In the numerical computations of flows around complex configurations, accurate 
calculations of force and moment coefficients for aerodynamic surfaces are required. 
When overset grid methods are used, the surfaces on which force and moment 
coefficients are sought typically consist of a collection of overlapping surface grids. 
Direct integration of flow quantities on the overlapping grids would result in the 
overlapped regions being counted more than once. The FOMOCO Utilities is a 
software package for computing flow coefficients (force, moment, and mass flow 
rate) on a collection of overset surfaces with accurate accounting of the overlapped 
zones. 

FOMOCO Utilities can be used in stand-alone mode or in conjunction with the 
Chimera overset grid compressible Navier-Stokes flow solver OVERFLOW. The 
software package consists of two modules corresponding to a two-step procedure: 
(1) hybrid surface grid generation ( MIXSUR module), and (2) flow quantities inte- 
gration ( OVERINT module). Instructions on how to use this software package are 
described in this user’s manual. Equations used in the flow coefficients calculation 
are given in Appendix A. 

The authors would like to thank the following for comments and discussions on 
the development of FOMOCO Utilities and for suggestions on the contents of this 
manuscript: Ray Gomez and Jim Greathouse from NASA Johnson Space Center, 
Steve Krist from NASA Langley, John Wai from Boeing, Lie-Mine Gea and Jeff 
Slotnick from McDonnell Douglas, Karen Gundy-Burlet, Earl Duque and our other 
colleagues at NASA Ames Research Center. 


* MCAT, Inc., San Jose, California 



1 Introduction and Overview 

The numerical computations of flows around complex configurations require ac- 
curate calculations of force and moment coefficients for aerodynamic surfaces. When 
overset grid methods (ref. 1) are used, the surfaces on which force and moment co- 
efficients are sought typically consist of a collection of overlapping surface grids (see 
figure 1 for example). Direct integration of flow quantities such as pressure and vis- 
cous stress on the overlapping grids would result in the overlapping regions being 
counted more than once. The FOMOCO Utilities is a software package for com- 
puting flow coefficients (force, moment, mass flow rate) on a collection of overset 
surfaces with accurate accounting of the overlapping zones. 

In this user's manual, instructions on how to run FOMOCO Utilities are given 
from a user’s point of view. Theoretical development of the methods used are 
described in detail in reference 2. The equations used for computing the flow coef- 
ficients are given in Appendix A. 



Figure 1. Overset surface grids with domain connectivity iblanks for wing/body /collar 
configuration (only every 3 points shown in streamwise direction for wing and 
collar). 


2 


The FOMOCO Utilities consist of two modules: the MIXSUR module and 
the OVERINT module. Both modules are written in Fortran 77. An overview of 
these modules is given in figure 2. The computation of flow coefficients on a set of 
overlapping surfaces is performed as follows. 

(1) A hybrid composite surface grid consisting of non-overlapping quadrilaterals 
and triangles is generated using the MIXSUR module. This is executed by typing 
mixsur < input parameters filename > output information filename 



Figure 2. Overview of FOMOCO Utilities and its modules. 
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(2) Flow coefficients are computed by integrating flow quantities on the hybrid 
surface using the OVERINT module. OVERINT can be used in stand-alone 
mode or in conjunction with the Chimera overset grid compressible Navier-Stokes 
flow solver OVERFLOW (ref. 3). In stand-alone mode, OVERINT is executed by 
typing 

overint < input parameters filename > output information filename 
In OVERFLOW mode, OVERINT is performed inside OVERFLOW. By setting 
the global parameter NFOMO to a positive integer in the OVERFLOW input 
parameters file, the flow coefficients are computed and reported every NFOMO 
steps. 

The files required to run FOMOCO Utilities are described in the subsequent 
sections. These include an input parameters file and the grid and solution files. 
Users who are not familiar with the PLOT3D (ref. 4) or OVERFLOW grid/solution 
file formats can refer to Appendix B for details. The same input parameters file is 
used to run both MIXSUR and OVERINT in stand-alone mode; and is described 
first in section 2. The MIXSUR and OVERINT modules are described in section 3 
and section 4, respectively. Concluding remarks are given in section 5. 


2 The Input Parameters File 

The same input parameters file is used for both MIXSUR and OVERINT in 
stand-alone mode. These parameters are used to 

(1) prescribe the flow conditions and the reference conditions used for flow coefficient 
computations, and 

(2) provide a description of the surface domain on which flow coefficient computa- 
tions are to be performed. 

A sample input parameters file is given in table 1(a) and brief descriptions of 
the parameters are given in table 1(b). More detailed discussions on the input 
parameters are provided in sections 2.1 and 2.2. 


2.1 Specification of Flow and Reference Conditions 

The flow conditions and reference conditions are specified in the first part of the 
input parameters file. The user should be careful in selecting the ALPHA (angle of 
attack) and BETA (side slip angle) parameters which are only used to transform the 
Cartesian force coefficients into lift, drag, and side force coefficients (see eq. (4.3) 
in section 4.3). 

For typical computational fluid dynamics (CFD) calculations, body axes are 
chosen such that the X-axis points nose to tail, the Y-axis points root to tip on the 
right wing, and the Z-axis points up. Let V^ and (uoo, v^. Wr^) be the free stream 
velocity magnitude and Cartesian components of free stream velocity, respectively. 
Then the angle of attack a and the side slip angle ft are defined as 

a=tair ‘fe)' (2i) 
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In cases where the body axes are aligned with the Cartesian axes as described 
above, the definition of a and 0 in eq. (2.1) results in the conventional lift, drag 
and side force coefficients. If the body axes are not aligned with the Cartesian axes, 
then the user must interpret the lift, drag and side force coefficients appropriately 
based on the ALPHA and BETA used in the input parameters file. In any case, 
the user can transform the Cartesian coefficients into coefficients under any other 
coordinate system. 

Different sets of reference conditions can be specified for each integration surface 
and component defined by the user (see section 2.2 for explanation of integration 
sin-face and component). The main point to note is that the reference condition 
parameters REFL, REFA, XMC, YMC, and ZMC should be in the same units as 
the input grid. 

Table 1(a). Sample input parameters file. 


0.8, 2.0, 0.0, 1.67e5, 1.4, 576.0 

FSMACH, ALPHA, BETA, REY, GAMINF, TINF 

2 

NREF 


1.0,871.8, 55.0, 0.0, 0.0 

REFL, REFA, XMC, YMC, ZMC } * 


1.0,168.4, 55.0, 9.0, 0.0 

}* 


3 

NSURF 


3,1 

NSUB, IREFS 

— 


1,3, 1,-1, 2, -2,1,1 

NG, IBDIR, JS, JE, KS, KE, LS, LE 



3, 3, 1,-1, 1,14, 1,1 




4, 2, 4, -1,1,1, 1,-1 



& 

2 

NPRI 



2,1 

NU1, NU2 } # 



3,1 

}# 

— 


2,2 

NSUB, IREFS 



2, 3,19,101,1,-1,1,1 

NG, IBDIR, JS, JE, KS, KE, LS, LE 


& 

3, 3,19,101,14, -1,1,1 




0 

NPRI 

— 


2,2 

NSUB, IREFS 



2, 3,101,-19,1,-1,1,1 

NG, IBDIR, JS, JE, KS, KE, LS, LE 


& 

3, 3,101,-19,14, -1,1,1 




0 

NPRI 

— 


2 

NCOMP 


WING 

“1 


2,2 

NIS, IREFC @ 


2,3 

ISC(N) % _| 


TOTAL 

“1 


3,1 

NIS, IREFC @ 


1,2,3 

ISC(N) % _| 



* Number of lines specified by NREF 

& Number of blocks specified by NSURF 

# Number of lines specified by NPRI 

@ Number of blocks specified by NCOMP 
% Provide NIS surface numbers 
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Table 1(b). Description of input parameters. The name of the solution file is q.save 
(see section 4.1). Parameters indicated by * should be in same units as input grid. 


FSMACH 


> 0 Free stream Mach number 



< 0 

Read free stream Mach number from q.save 

ALPHA in range [-360,360] 

Angle of attack in degrees 


outside of [-360,360] 

Read angle of attack from q.save 

BETA 

in range [-360,360] 

Side slip angle in degrees 


outside of [-360,360] 

Read side slip angle from q.save 

REY 

> 0 

Reynolds number per grid unit for viscous flow 


= 0 

Inviscid flow assumed 


< 0 

Read Reynolds number from q.save 

GAMINF > 0 

Free stream ratio of specific heats 


< o 

Read free stream ratio of specific heats from q.save 

TINF 

> 0 

Free stream temperature in Rankine 


< 0 

Read free stream temperature in Rankine from q.save 

NREF 


Number of sets of reference quantities 

REFL 

7^0 

Reference length for flow coefficients * 


= 0 

Use 1.0 as default reference length 

REFA 

^0 

Reference area for flow coefficients * 


= 0 

Use total integrated surface area from a component 


XMC,YMC,ZMC 

NSURF 

NSUB 

IREFS 

NG 

IBDIR 

JS,JE,KS,KE,LS,LE 


in which surface belongs, or if surface does not 
belong to any component, use area of surface 
X,Y,Z coordinates of moment axes center * 

Number of surfaces to compute coefficients 

Number of subsets that belong to the surface 
Reference quantity set number for surface 

Grid number 

Direction of' surface normal 

(1 = +J, -1 = — J, 2 = +K, -2 = -K, 3 = +L, -3 = -L) 
Start and end indices in J, K, and L, respectively. If an 
input index is negative, (IMAX+1) is added to form a 
positive index where IMAX is the grid dimension in I 
(I = J, K, or L), e.g., JS = —1 is the same as JS = JMAX. 


NPRI 

NU1,NU2 

NCOMP 

CNAME 

NIS 

IREFS 

ISC(N) 


Number of subset pairs to specify priority 

Subset numbers within surface where subset NU1 will 

be kept if it overlaps subset NU2 

Number of components (can be omitted if no component) 
Component name 

Number of integration surfaces for this component 
Reference quantity set number for component 
List of surface numbers of integration surfaces 
for this component 
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2.2 Specification of Surface Domain 

The surface domain on which flow coefficients are computed is specified in the 
input parameters file. A complex surface domain is typically divided into one or 
more ‘integration surfaces’ where each integration surface consists of a collection 
of subsets from different grids in the input grid file grid. in (see section 3). Each 
subset can be any part of a J, K, or L=constant slice in the computational domain 
of any grid. While the flow solution is assumed to be node based, the integration 
scheme is cell (face) centered. For this reason, the indices used to specify the subsets 
(JS,JE,KS,KE,LS,LE) should include the end points of the region to be integrated. 
Grids in the input grid file that do not contain any of the specified subsets are 
disregarded. 

It is frequently convenient to group a collection of integration surfaces together 
to form a ‘component’. Flow coefficients for each component and each integration 
surface are computed by the integration module OVERINT. Each component and 
each integration surface is allowed to have a different set of reference quantities 
(length, area and moment center). The hierachy of surface domain, components, 
integration surfaces and subsets is shown in figure 3. Note that an integration 
surface may belong to more than one component although this is not typically the 
case. 



| j = Surface Subset 

Figure 3. Hierachy of surface domain, components, integration surfaces, and sub- 
sets. 


For example, a wing/body surface domain can be divided into an upper wing 
integration surface, a lower wing integration surface and a body integration surface 
where each integration surface consists of several subsets from different grids. The 
upper wing and lower wing may be grouped into a ‘wing component’. 
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In the overlapped zone between two subsets from the same integration surface, 
the code automatically determines which subset has fewer points in the overlapped 
zone. Points from this subset are then blanked in the overlapped zone. The user 
has the option to override this scheme by specifying priorities for pairs of subsets. 
This is achieved by making NPRI > 0, and then providing two subset numbers for 
each pair (see table 1(a)). Points from the second subset will be blanked if they are 
in the overlapped zone with the first subset (see section 3.3 for an example). 

Triangles are automatically constructed in the gap produced between neighbor- 
ing subsets after the blanking step. These triangles are typically one to two layers 
deep (zipper grids). The blanking of points and construction of triangles are not 
considered between subsets belonging to different integration surfaces even if the 
subsets overlapped or abutted. Hence, if the user does not wish to have triangles 
created between two abutting subsets, the two abutting subsets should be defined 
under two different integration surfaces. 

The surface subsets specified in the input parameters file may belong to solid 
surfaces (no flow-through) or field surfaces (flow-through). For example, the surface 
domain for a nacelle calculation may include the nacelle walls (solid) and the inlet 
and outlet surfaces (field). For a solid surface, the momentum flux and mass flow 
rate through the surface are zero. 


3. The MIXSUR Module 

The MIXSUR module is used to create a hybrid composite grid consisting of 
non-overlapping quadrilaterals and triangles on which flow quantities are integrated. 
The hybrid grid is automatically generated from the information contained in the 
input files described in section 3.1 below. This procedure can be performed prior to 
computing the flow solution since only the volume grids with domain connectivity 
iblanks are required as input. 

The user is urged to use the automatically generated PLOT3D (ref. 4) command 
files to view and check the integrity of the composite hybrid grid before proceeding 
to the next step (see section 3.2). Note that a comparable grid resolution in the 
overlapped region between grids can greatly enhance the success of the triangulation 
of the region. Moreover, this property is highly desirable for accurate information 
transfer between grids. 


3.1 Input and Output Files 

The MIXSUR module requires two input files. The first is an input parameters 
file (see section 2) which describes the surface domain on which flow coefficients are 
to be computed. The second is a PLOT3D single or multiple volume grid file with 
domain connectivity iblanks. This file should be named grid. in and is identical to 
the input file required by the flow solver OVERFLOW. The user should be aware 
of the following two points on the grid. in file. 

(1) For viscous computations, an alternative to supplying the entire volume grid for 
each grid is to provide a volume grid with two surface slices: the grid surface 
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on which flow coefficients are to be computed and the neighboring grid surface 
in the normal direction. If only inviscid flow is considered, an alternative to 
supplying the volume grids is to provide just the grid surfaces on which flow 
coefficients are to be computed. 

(2) The domain connectivity iblanks in grid. in are only used to determine the 
blanked cells on the surface domain. The MIXSUR code has its own methods 
to compute grid connectivity and does not rely on the connectivity information 
contained in the iblanks in grid. in. If grid. in does not contain iblanks, MIXSUR 
will assume a value of 1 for all the domain connectivity iblanks. 

Several output files are written by MIXSUR. Some of these files are used as 
input to the OVERINT module while others are used for visualization of the hybrid 
grid produced by MIXSUR. Summaries of the input and output files are given in 
tables 2(a) and 2(b). 


Table 2(a). Input files required by MIXSUR. 


Filename 

Description 

STDIN 

Input parameters in ASCII format. 

Grids not specified in this file but in grid. in are disregarded 
by MIXSUR and OVERINT. 

grid, in 

Input grids in unformatted PLOT3D volume grid file with domain 
connectivity iblanks (single or multiple grid format). 


Table 2(b). Output files produced by MIXSUR. 


Filename 

Description 

mixsur.fmp 

ASCII file containing parameters of MIXSUR that will be read by 
the flow solver OVERFLOW for running OVERINT. 

grid.ibi 

Surface subsets in unformatted PLOT3D grid file with integration 
iblanks (single or multiple grid format). Each grid is a surface grid 
defined by the J, K or L surface subsets in the input parameters 
file. The integration iblanks indicate which cells are to be integrated. 

grid, ns f 

Surface subsets in unformatted PLOT3D grid file (single or multiple 
grid format). Each grid is a surface grid next to the corresponding 
grids in grid.ibi in the user-defined normal direction. This file is 
generated if grid. in contains more than a single slice in the first grid. 

grid.bnd 

Unformatted PLOT3D multiple grid file. Each grid is a 1-D string of 
points marking the boundary of the region to be triangulated. 

grid.tr i 

Unformatted PLOT3D multiple grid file. Each grid is a 1-D string of 
points marking the corners of each triangle in the triangulated region.* 

grid. nor 

Unformatted PLOT3D multiple grid fake solution file. The 2nd, 3rd 
and 4th component of each solution contains the x,y,z components of 
the normals of the triangular cells in grid.tri.* 

grid.ptv 

File containing pointers to each vertex of each triangle.* 


* Generated only if triangles are needed. 
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3.2 Grid Visualization 

The hybrid composite grid can be viewed using PLOT3D on the files grid.ibi 
and grid.tri. Table 3 shows the PLOT3D command files that are automatically 
written out. Some example plots using these command files are given in section 3.3. 


Table 3. PLOT3D command files generated by MIXSUR. 


Filename 

Description 

ib.com 

Surface subsets with integration iblanks (quadrilaterals) 

bnds.com 

Gap boundary strings 

zips.com 

Zipper grids (triangles) 

zn.com 

Normals for zipper grids 

nowalls.com 

Remove all walls from any of above 

minmax.com 

Bounding box for all grids 


To view the hybrid composite grid followed by normals of the zipper grids: 

(1) Type @ib at the PLOT3D prompt. Type return after the plot is up. 

(2) Type @zips at the PLOT3D prompt. Type return after the plot is up. 

(3) Type @zn at the PLOT3D prompt. 

To view the boundary strings of the triangulated regions: 

(1) Type @ib at the PLOT3D prompt. Type return after the plot is up. 

(2) Type @bnds at the PLOT3D prompt. 


3.3 Examples 

A simple configuration consisting of three grids is considered: the body, the wing 
and the collar (ref. 5). The overlapping surface grids on which flow coefficients are 
computed are shown in figure 1. Figure 4(a) shows the overlapping surface grids 
on the body/collar junction that are used as input to MIXSUR. Note that part 
of the body grid has been blanked out due to the presence of the wing. This 
blanking information is typically generated by a domain connectivity program such 
as PEGSUS (ref. 6) or DCF3D (ref. 7) and is stored in the grid. in file. 

Figure 4(b) shows the surface grids after the integration iblank array has been 
constructed by MIXSUR. Gaps are created between neighboring grids while the 
remaining unblanked points belong to quadrilaterals that do not overlap. This view 
can be generated by using the ib.com command file in PLOT3D. 

Points along the boundaries of the gaps are identified and connected to form gap 
boundary strings. These strings can be viewed by using the bnds.com command 
file in PLOT3D as shown in figure 4(c). 
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(b) 

Figure 4. Body /collar junction, (a) Overlapping surface grids with domain connec- 
tivity iblanks, (b) surface grids with integration iblanks. 
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(d) 

Figure 4. Body/collar junction, (c) Gap boundary strings, (d) hybrid composite 
surface grid with non-overlapping quadrilaterals and triangles. Collar has higher 
priority than body (automatically determined by default). 
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Figure 4. Body/collar junction, (e) Hybrid composite surface grid where the body 
subset is manually given higher priority over the collar subset. 

Finally, triangles are constructed to fill the gaps. The resulting hybrid grid can 
be viewed by using ib.com followed by zips.com in PLOT3D as shown in figure 4(d). 
With this view up in PLOT3D , the normals of the triangular cells created can be 
checked by using the zn.com command file. Since the MIXSUR code is known to 
produce some bad triangles in difficult cases, the user is strongly urged to use these 
PLOT3D command files to check the integrity of the hybrid grid created prior to 
running the OVERINT module. 

The hybrid grid in figure 4(d) is generated based on the default blanking priority, 
i.e. points from the coarser of two overlapping subsets are blanked. In this case, 
points from the body are blanked in the overlapped zone by the default scheme. 
The default priority may be manually overidden by the user for reasons discussed 
in paragraph (c) in section 3.4. In this example, the priority of the body subset can 
be made to be higher than that of the collar subset. The resulting hybrid grid is 
shown in figure 4(e). 

An example of a complex configuration is demonstrated by the External Tank 
(ET) of the Space Shuttle Launch Vehicle as shown in figure 5. The view is looking 
downstream from the interior of the ET base region. The ET surface grid is over- 
lapped by the collar surface grids from the various feedlines that intersect the ET. 
Further examples of complex configurations can be found in reference 2. 


Figure 5. Hybrid composite grid for the base region of the External Tank. 


3.4 Special Notes on Using MIXSUR 

Given the input grids and the specification of the overlapping surface subsets, 
the hybrid composite grid is automatically generated by MIXSUR. High quality 
triangles are usually created when the grid resolutions of neighboring grids are 
comparable in the overlapped zones, and when the surface geometry is sufficiently 
resolved by the grid. However, when there is a large discrepancy in grid resolution 
between neighboring grids in the overlapped zones, the triangles created are not 
guaranteed to be totally correct. A triangle is considered ‘bad' if it overlaps an 
unblanked quadrilateral. Its surface normal is pointed in the opposite direction to 
neighboring good triangles. 

Figures 6(a) and 6(b) show the hybrid grid between two overlapping subsets and 
the surface normals of each triangle. For each triangle, the normals are plotted at 
the vertices (i.e., three arrows pointing in the same direction are plotted for each 
triangle). A bad triangle is produced in the region of high grid resolution mismatch. 
Its surface normal is pointing in the opposite direction to its good neighbors. 


Also shown in figure 6 is a singular axis point in one of the subsets. Such 
a point is surrounded by triangular cells and may be found at a domain bound- 
ary in a structured quadrilateral mesh. Singular axis points are automatically de- 
tected by MIXSUR and the neighboring triangular cells are treated as a zipper grid, 
i.e., OVERINT will use the integration scheme for triangles on these cells. 



(a) (b) 

Figure 6. Hybrid composite grid for two overlapping subsets. A bad triangle indi- 
cated by thick lines is shown in the region of high grid resolution mismatch where 
the normals are pointing in the wrong direction, (a) Top view, (b) isometric view 
with surface normals of triangles (zipper grid). 

The following notes may help to reduce the number of bad cells created: 

(a) A large number of surface subsets describing a complex configuration should 
be split into multiple integration surfaces where subsets belonging to different 
integration surfaces do not communicate with each other. 

(b) High curvature regions or extremely sharp corners are sometimes trouble regions. 
For example, points on the lower surface of a wing may be connected to the 
upper surface near the trailing and/or leading edge. This problem can be fixed 
by defining two integration surfaces for the wing: an upper surface with its 
collection of surface subsets and a lower surface with its collection of surface 
subsets. Now, points on the upper and lower surfaces are not allowed to be 
connected to each other. 

(c) The use of the default blanking priority between overlapping subsets (points 
from the coarser subset is blanked) can usually result in good quality triangles 
in the hybrid grid. However, this is not always true in cases where there is severe 
grid resolution mismatch in the overlapped zones. For these cases, a smoother 
triangulation gap boundary may sometimes be obtained by manually imposing 
blanking priorities using the NPRI, NU1, and NU2 parameters (see section 2). 
This may then result in better quality triangles. For example, the hybrid grid for 
the body /collar junction in figure 4d was generated using the default priorities 
which produced smooth gap boundaries for the triangulation. However, if the 
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body is given priority over the collar, a more rugged gap boundary is generated 
as shown in figure 4(e). In this case, the triangles generated are still of good 
quality, but the results from using the default priorities are preferred since the 
gap boundaries are smoother. 

(d) The blanking scheme assumes that surface grids in the overlapped zones are 
approximating the same surface geometry, i.e., points from any one grid in the 
overlapped zone should lie very close to the bilinear surface defined by the net- 
work of points from any other grid in the overlapped zone. If ‘big’ gaps exist 
between grids in the overlapped zone, some points may not be blanked and the 
subsequent triangulation will be bad. 


4. The OVERINT Module 

The OVERINT module is used to compute and report the force, moment, and 
mass flow rate coefficients on the hybrid composite grid generated by the MIXSUR 
module. There is no interpolation of the flow solution onto new grid points since 
none are created by MIXSUR. 

OVERINT can be run in stand-alone mode or be called from the flow solver 
OVERFLOW. In stand-alone mode, flow coefficients at the time step level for the 
given solution file are reported. In OVERFLOW mode, a history of the flow coeffi- 
cients versus time step and versus elapsed CPU time is reported in an OVERFLOW 
output file. Software utilities for extracting and plotting histories of particular co- 
efficients are available under the OVERFLOW release package. 

When using OVERINT in OVERFLOW mode for moving body problems, 
MIXSUR still only needs to be run once before the flow solver, provided the grids 
do not deform. This is because grid points on the surface grids of rigid moving com- 
ponents remain in constant relative positions to each other. Hence, the connections 
between grid points on the hybrid surface grid and the integration iblank informa- 
tion produced by MIXSUR remain unchanged during body motions, although the 
physical coordinates of the grid points are changing. The connection and blanking 
information from MIXSUR is read by OVERINT which then uses the most recent 
grid coordinates to compute cell areas and normals at each time step. 


4.1 Input and Output Files 

Five input files are required to run OVERINT in stand-alone mode (see ta- 
ble 4(a)). The first is the same input parameters file that is used to run MIXSUR. 
The next three files, grid.ibi, grid.nsf, and grid.ptv, are files produced by MIXSUR. 
The final file is the solution file on the grid system given in grid.in. This file is as- 
sumed to be named q.save and can be in OVERFLOW or PLOT3D solution file for- 
mat. OVERINT automatically detects whether the solution file is in OVERFLOW 
or PLOT3D format and reads it appropriately. Note that an OVERFLOW solu- 
tion file may contain more information than a standard PLOT3D solution file. The 
additional information may include parameters for species convection, dependent 
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variables for turbulence tranport equations and species convection equations (see 
Appendix B for details of the formats). 

Three input files are required to run OVERINT in OVERFLOW mode (see 
table 4(a)). The grid.ibi and grid.ptv files are the same as those used by the code 
in stand-alone mode. The third file needed is the mixsur.fmp file generated by the 
MIXSUR module. 

Table 4(a). Input files required by OVERINT. 


Filename 

Description 

STDIN 

Input parameters file in ASCII format 
(same as MIXSUR ’s input) (s) 

mixsur.fmp 

(generated by MIXSUR ) (o) 

grid.ibi 

(generated by MIXSUR ) ( s,o ) 

grid, ns f 

(generated by MIXSUR ) (s) 

grid.ptv 

(generated by MIXSUR ) (s,o) 

q.save 

OVERFLOW or PLOT3D solution file * (s) 


* (single or multiple grid format) 

(s) Required by stand-alone mode. 
(o) Required by OVERFLOW mode. 


A different output file is produced by OVERINT in stand-alone mode and 
OVERFLOW mode (see table 4(b)). 

Table 4(b). Output file produced by OVERINT. 


Mode 

Filename 

Description 

Stand-alone 

STDOUT 

flow coefficients, integrated areas 

OVERFLOW 

fomoco.out 

history of flow coefficients, integrated areas 


4.2 Input Data 

The dependent flow variables from the q.save solution file are assumed to be 
non-dimensionalized in the same manner as in the OVERFLOW flow solver (see 
eq. (A. 3. 2b) in Appendix A). If a spatially variable ratio of specific heats is given in 
q.save, that value will be used. Otherwise, the free stream value GAMINF in the 
input parameters file is used (see table 1). 

The following parameters in the input parameters file, REFL (reference length), 
REFA (reference area), XMC, YMC, ZMC (moment reference center), should be 
provided in the same units as the input grid in grid.in. Similarly, the input Reynolds 
number REY is the Reynolds number based on unit length of the grid system. For 
example, if the grid system is represented in feet, REY should be the Reynolds 
number per foot. 
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4.3 Output Data 

The OVERINT code computes projected areas, force, moment, and mass flow 
rate coefficients for each integration surface and component defined. Each integra- 
tion surface and component can have its own set of reference quantities (length, area, 
moment center) from which the coefficients are computed. In stand-alone mode, co- 
efficients for each integration surface and component are reported. In OVERFLOW 
mode, only the coefficients for each component are reported. Table 5(a) shows a 
sample output data file for OVERINT in stand-alone mode. Table 5(b) shows the 
pseudo code used to write the flow coefficient histories file when using OVERINT 
in OVERFLOW mode. The data reported include: 

(1) Projected areas in the X, Y, and Z directions and the total area. 

(2) X, Y, and Z force coefficients (pressure, viscous, momentum flux, and total). 

(3) Lift, drag, and side force coefficients (pressure, viscous, momentum flux, and 
total). 

(4) Moment coefficients about X (roll), Y (pitch), and Z (yaw) axes centered at 
(XMC, YMC, ZMC). 

(5) Mass flow rate coefficient. 

In stand-alone mode, the momentum flux and mass flow rate coefficients are 
reported only if they are non-zero; and X, Y, Z force coefficients are reported for 
components only. 

Equations used to compute the non-dimensional force F, moment A4 and mass 
flow rate m are given in Appendix A. The force coefficients Cf, moment coefficients 
C m and mass flow rate coefficient C p are obtained as 


C f = 


T 

Qoo A-ref 


c m = 


M 

Q oo^ief Aref 


c a 


m 


Poo Lqo A r ef 


(4.1) 


where 

Qo o = 2 Poo Loo = 2 P°° M OC ^oo (4-2) 

and (') is used to denote non-dimensional variables (see section A3 in Appendix A); 
the subscript <*, is used to denote free stream quantities; p, V, M, and c are the 
density, flow speed, Mach number and sound speed, respectively; Z re f and A re f are 
the reference length (REFL) and area (REFA), respectively. 

The lift, drag and side forces ( Ci,Cd,C s ) are computed from the X, Y, and Z 
components of forces (C x , C y . C z ) as follows. 


Ci = — C x sin a + C z cos a 

Cd — {C x cos a + C Z sin a) cos (3 — C y sin / 3 (4.3) 

C s = ( C x cos a + C z sin a) sin (3 + C y cos / 3 


The roll, pitch, and yaw moment coefficients are counted positive for right hand 
rotations about the X, Y, and Z (body) axes, respectively. 
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Table 5(a). Sample output data file from OVERINT in stand-alone mode 


Output from OVERINT version 

1.0 

Mach number * 
Beta = 
Gamma_inf = 

0.800 Alpha = 2.000 

0-000 Reynolds number * 
1.400 Temperature„inf * 

0. 16700E+06 
76.000 

Reference Quantities set number 1 

Reference length for coefs * 0.10000E+01 

Reference area for coefs * O.B7180E+03 

Moment axes center=0 . 55000E+02 0.00000E+00 O.OOOOOE+OO 

Reference quantities set number 2 

Reference length for coefs = 0.10000E+Q1 

Reference area for coefs * 0.16840E+03 

Moment axes center*0 . 55000E+02 0.90000E+01 O.OOOOOE+OO 

Data for surface 1 



Integrated Areas 

x y 

Z 

Total 

-0.352421X4-01 

0 . 444506E+03 

0 . 162479E 

-04 0 . 7 03 312E+03 

Force coefs 

Lift 

Drag 

Side 

Pressure 

Viscous 

Total 

0.005121 

-0.000037 

0.005085 

0.001412 

0.002498 

0.003911 

0.029349 

0.000040 

0.029389 

Moment coefs 

Roll (X) 

Pitch (Y) 

Yaw (Z) 


0.000131 

0.170961 

0.264587 

Data for surface 2 



Integrated Areas 

x y 

z 

Total 

0 . 222120E-06 

0.1849162+01 

-0 .828426E+02 0.842213E+02 

Force coefs 

Lift 

Drag 

Side 

Pressure 

Viscous 

Total 

-0.031813 

-0.000139 

-0.031952 

0.003530 

0.002181 

0.005711 

-0.002711 

0.000018 

-0.002693 

Moment coefs 

Roll <X) 

Pitch <Y) 

Yaw (Z) 


-0.008758 

■0.825687 

0.078427 

Data for surface 3 



Integrated Areas 

X Y 

z 

Total 

0 . 2634 19E-06 

0.184916E+01 0 . 828425E+02 0.842213E+02 

Force coefs 

Lift 

Drag 

Side 

Pressure 

Viscous 

Total 

0.137826 

0.000103 

0.137929 

0.001455 

0.002138 

0.000683 

0.006899 

-0.000038 

0.006861 

Moment coefs 

Roll (X) 

Pitch (Y) 

Yaw (Z) 


0.082796 

3.609517 

-0.166994 


Data for component WING 

Integrated Areas 

X Y 

z 

Total 

0 . 48 5539E-06 

0.369832E+01 

. -0 . 15258 BE 

-04 0.168443E+03 

Force coefs 

X 

Y 

Z 

Pressure 

-0.001626 

0.004188 

0.106021 

viscous 

0.004318 

-0.000020 

0.000115 

Total 

0.002691 

0.004168 

0.106136 

Force coefs 

Lift 

Drag 

Side 

Pressure 

0.106013 

0.002075 

0.004188 

Viscous 

-0.000035 

0.004319 

-0.000020 

Total 

0.105978 

0.006394 

0.004168 

Moment coefs 

Roll (X) 

Pitch (Y) 

Yaw (Z) 


0.074037 

2.783830 

-0.088568 

Data for component TOTAL 



Integrated Areas 



X 

Y 

Z 

Total 

-0 . 352421E+01 

0.448204E+03 

1 O.OOOOOOE+OO 0 . 87 1755E+03 

Force coefs 

X 

Y 

Z 

Pressure 

0.000919 

0.030158 

0.025647 

Viscous 

0.003332 

0.000036 

0.000073 

Total 

0.004251 

0.030194 

0.025720 

Force coefs 

Lift 

Drag 

Side 

Pressure 

0.025599 

0.001813 

0.030158 

viscous 

-0.000044 

0.003333 

0.000036 

Total 

0.025556 

0.005146 

0.030194 

Moment coefs 

Roll (X) 

Pitch ( Y) 

Yaw (Z) 


0.198947 

0.708695 

0.242800 
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Table 5(b). Pseudo code used to write flow coefficient histories when using OVERINT 
in OVERFLOW mode. 


DO 1=1 .NITER 

IF ( (MOD(I,NFOMO)=0) OR (l=NITER) ) THEN 
DO NC=1,NCOMP 
WRITE CNAME(NC) 

WRITE l,-NC,SAC(NC),TAXC(NC),TAYC(NC),TAZC(NC) 

WRITE CXPC(NC),CYPC(NC),CZPC(NC),CXVC(NC),CYVC(NC),CZVC(NC) 
WRITE CXMC(NC),CYMC(NC),CZMC(NC),CLPC(NC),CDPC(NC),CSPC(NC) 
WRITE CLVC(NC),CDVC(NC),CSVC(NC),CLMC(NC),CDMC(NC),CSMC(NC) 
WRITE CMRC(NC),CMPC(NC),CMYC(NC),CFRC(NC),TIMECP(I) 

END DO 
END IF 
END DO 


CNAME 

NITER 

NFOMO 

I 

NC 

S AC,T AXC.T AY C.T AZC 

CXPC.CYPC.CZPC 

CXVC.CYVC.CZVC 

CXMC.CYMC.CZMC 

CLPC.CDPC.CSPC 

CLVC.CDVC.CSVC 

CLMC.CDMC.CSMC 

CMRC.CMPC.CMYC 

CFRC 


: component name 
: number of time steps in flow solver 
write force/moment data every NFOMO steps 

: time step number 
: component index number 

total surface area, projected X, Y, Z areas 
X, Y, Z pressure force coefficients 
X, Y, Z viscous force coefficients 
X, Y, Z momentum flux coefficients 
lift, drag, side pressure force coefficients 
lift, drag, side viscous force coefficients 
lift, drag, side momentum flux coefficients 
roll, pitch, yaw moment coefficients 
mass flow rate coefficient 


TIMECP 


= elapsed CPU time 


4.4 Special Notes on Using OVERINT 

(1) When computing flow coefficients for a single grid, the MIXSUR module still 
needs to be run before the OVERINT module. 

(2) The OVERINT code integrates (p — poo) over the defined subsets where p is the 
pressure. If the defined subsets form a closed surface, then the true aerodynamic 
force and moment coefficients are computed. If the defined subsets form an open 
surface, the user must decide whether to adjust the reported coefficients by the 
poo term. 

(3) The user should be aware that the answers obtained with the current flow quan- 
tity integration scheme (see Appendix A) may vary slightly depending on the 
following: 

(a) Quadrature rule used to integrate a function over the quadrilaterals and 
triangles. 

(b) Integration using triangles versus quadrilaterals. 

(c) Size and shape of the triangulated region. 

Initial results seem to indicate that the variations caused by the above factors 
are typically much smaller than the accuracy of the flow coefficients required. 
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5 Concluding Remarks 

The FOMOCO Utilities are currently in beta release stage. The software has 
been used by a number of groups within NASA and industry. These include oblique 
wing with pylons, fins and nacelles (ref. 2), subsonic transport with wing, body, 
pylon, nacelle, core-cowl (refs. 2 and 8), Stratospheric Observatory for Infrared 
Astronomy (SOFIA) (ref. 9), Space Shuttle Launch vehicle (refs. 2 and 10), and 
the X-CRV crew return vehicle (ref. 10). Very good results have been obtained for 
comparisons between numerical computations and experiments. 

Bug reports and comments are most appreciated and can be sent to 
wchan@nas.nasa.gov 
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Appendix A 

The equations used to compute forces, moments, and mass flow rate for a surface 
are described in this appendix. 

A1 Surface Integrals of Flow Quantities 

In general, the forces and moments acting on a solid or field surface consist of 
contributions from the pressure, viscous stress and momentum flux. For a solid 
surface, the momentum flux and mass flow rate through the surface are zero. 

In Cartesian tensor notation, the force Ti acting on a surface A due to a moving 
fluid is given by 

J- t = I aijfij dA (-4.1.1) 

J A 

where <r jj is the stress tensor and n* is the local unit normal to the surface. The 
force 0i due to momentum flux through a field surface is given by 

0; = / pUiUjfij dA (A. 1.2) 

Ja 

where p is the density and u, is the velocity. The mass flow rate m through a field 
surface is given by 

m = / pujTij dA (A. 1.3) 

Ja 


A2 General Form of the Stress Tensor 

The stress tensor is given by 

C ij — p&ij ”1” Tij (A. 2.1) 


where p is the pressure, S lJ is the Kronecker delta function and r tJ is the viscous 
stress tensor given by 


Tij = n 




+ A Sij 


duk 

dx k 


(A. 2.2a) 


The laminar viscosity is denoted by p. From the Stokes hypothesis A = -§/q 
eq. (A. 2. 2a) becomes 


Tij = p 




2 duk 

3 tj dx k 


(A.2.26) 


When the stress tensor is evaluated at a solid wall, the turbulent eddy viscosity can 
be assumed to be zero. 
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A3 Non-dimensional Form of the Stress Tensor 

Let (~) denote non-dimensional variables and ( ) (X) denote the free stream value 
of the variable. A typical term in the stress tensor is given by 


du 

a p + p— 

dx 

The non-dimensionalization used in OVERFLOW is summarized by 

x d _ 1 d 
Z’ dx = Ld± 

P - P 


x = 


it 


(A. 3.1) 

(A. 3. 2a) 
(A. 3. 2b) 

Coo Poo poo 

where L is the length scale and c is the sound speed. Substituting in eq. (A. 3.1) 
gives 

// . . ft ri( r _ _ ft \ 

(A. 3.3a) 

±j UJb 

or 


u = 


P = 


P = 


P = 


P<x>Coc 


a — 


P OC C 00 


(T ~ — p + 


L dx 

Poo du 


poo Coo L dx 


(A.3.36) 


The Reynolds number Re and the scaled Reynolds number Re used internally in 
OVERFLOW are defined as follows 


Re = 


Poo Loo -L 

Poo 


Re = 


Poc Coo L Re 

Poo Afoo 


(A.3.4) 


where M 0 c is the free stream Mach number. Using the scaled Reynolds number in 
eq. (A.3.4) and dropping all hats, the non-dimensional form of the stress tensor is 
given by 

_ p ' 

&ij — pOij V — — 

Re . 


dui duj\ 2 duk 

dxj dxi ) 3 13 dxk 


(A. 3.5) 


A4 Stress Tensor in Cartesian and Generalized Coordinates 


In Cartesian coordinates, the components of the non-dimensional stress tensor 
(following from eq. (A. 3. 5)) are given by 


& XX 

0 r yy 
<7zz 
® xy 
<Xyz 


2 p 

-p + ~^{2u x - Vy - w z ) 

3Re 

2p . 

-p + —=r(2vy -w z - u x ) 

ZRe 

2 p 

~P + ^=r(2u> z - U x - Vy) 

3 Re 

P / \ 

\Uy + V x ) 

Re 

P 


Re 


{v z + Wy) 


®XZ — 


p / \ 

— {U z + W x ) 

Re 


(4-4.1) 
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The right hand sides of eq. (A.4.1) can be expressed in generalized coordinates 
77 , and (, using the chain rule. For example, a xx can be written as 


— [2(w^ x d" UfjTjj; + U^ x ) “b VyVy ”b 

3 Re 


G xx — P ”b 

— (uie£r -j- -j/i n -1- I/!,/ - II 


(A 441 


Quantities on the right hand side of eq. (A. 4. 2 ) can be conveniently computed in 
the flow solver. 


Typically, forces and moments acting on a solid wall are required. Suppose the 
solid wall is along a £ = constant surface. Since all velocity components are zero 
at a viscous solid wall, all terms containing derivatives of velocity components in £ 
and 77 are exactly zero in the viscous part of the stress tensor when it is evaluated at 
the wall. Note that no thin layer approximation is invoked in dropping the £ and 77 
derivative terms. The components of the non-dimensional stress tensor at the wall 


now become 

2 p . 

<?xx = ~P + — (Cz«c - m ) 

Re 

2 /z_ v 

Gyy = ~P+ 

Re 

2 /x . 

o zz = -p+ — - m) 

Re 

fJL 

&xy ~ H“ Cy^C/ 

Re 

°yz = + <>c) 

Re 

fi 

O'xz — ”b G«c) 

Re 


(A4.3a) 


where 

m = -(Cx«c + Cy«c + Cz™<) 


(A.4.36) 


Consider the terms involving derivatives of velocity in the tangential direction 
to the surface in eq. (A. 4. 2). There are two situations when these terms may be 
non-zero. The first is in the case of a field (flow-through) surface where tangential 
velocity gradients are not necessarily zero. The second is in the case of a solid 
surface in unsteady motion. In an inertial reference frame fixed in space, different 
points on the solid surface may have different velocities (e.g., a rotating body), 
thus resulting in non-zero tangential velocity gradients. In the above two cases, 
all derivative terms must be retained in computing the stress tensor components. 
These tangential derivative terms are disregarded in the current release of FOMOCO 
Utilities. They will be implemented in the near future. 
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A5 Flow Quantities Integration 

The total force Ti, moment Mi and mass flow rate m are computed by summing 
the respective elemental forces A Ti, elemental moments AMi and elemental mass 
flow rates Am. The elemental quantities are computed in a cell-centered scheme 
for each cell in the surface domain consisting of non-overlapping quadrilaterals and 
triangles. 

The elemental force A Ti acting on a vector area element tiiAA where n* is the 
local unit normal at the surface is given by 

A Ti — aijTijAA (A. 5.1) 

where (“) denotes some averaging operator (see below). Similarly the elemental 
moment A Mi about a point rf is given by 

A Mi = t ijk {r j - r°)AF fc (A.5.2) 

where r j is the position vector of the centroid of the local area element and tij k is 
the Levi-Civita permutation symbol. The elemental mass flow rate Am is given by 

Arrt = pTIjnjAA (A. 5. 3) 

The ( — ) quantities above are evaluated at the cell center using a zeroth order 
scheme where the value at the center is taken to be the arithmetic average of the val- 
ues at the vertices. Higher order schemes such as assuming a bilinear quadrilateral 
or triangle may be used but are not employed here. Simple numerical experiments 
seem to indicate that such higher order schemes are not needed for the typical 
accuracy required of the flow coefficients. 


A6 Evaluation of Area Elements and Metrics 

The local vector area element UiAA for a quadrilateral is computed by one half 
times the cross product of its diagonals. Similarly, the local vector area for a triangle 
is computed by one half times the cross product of the two vectors along any two 
of its edges. In both cases, the normal is ensured to be outward pointing from the 
surface. 

The metric terms (( x ,Cy,(z) — VC are evaluated by 


VC = rnAA/AV (A. 6.1) 

where AV is the local cell volume. This local cell volume is evaluated in a finite 
volume sense as in references 11 and 12. However, a symmetric evaluation indepen- 
dent of the index direction is used. The symmetric value is obtained by averaging 
the four volumes computed using the four diagonals of the hexahedral cell. 
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Appendix B 

The various file formats used by FOMOCO Utilities are given below. 


PLOT3D multiple grid XYZ file with IBLANKS 

READ(IUNIT) NGRID 

READ(IUNIT) (JDIM(N),KDIM(N),LDIM(N),N=1 , NGRID) 

DO 10 N=1, NGRID 

READ(IUNIT) (((X(J,K,L),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)), 

& (((Y(J,K,L),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)), 

& (((Z(J,K,L),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)), 

& (((IBLANK(J,K,L),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)) 

10 CONTINUE 


PLOT3D multiple grid solution file 

READ(IUNIT) NGRID 

READ(IUNIT) (JDIM(N),KDIM(N),LDIM(N),N=1 , NGRID) 

READ(IUNIT) FSMACH, ALPHA, REY, TIME 
DO 10 N=1, NGRID 

READ(IUNIT) ((((Q(J,K,L,NV),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)),NV=1 ,5) 
10 CONTINUE 


OVERFLOW multiple grid solution file (most general form since version 1 .7e) 


READ(IUNIT) NGRID 

READ(IUNIT) (JDIM(N),KDIM(N),LDIM(N),N=1 ,NGRID),NQ,NQC 
READ(IUNIT) FSMACH, ALPHA, REY, TIME, GAMINF, BETA, TINF, 

& IGAM,HTINF,HT1 ,HT2,(RGAS(I),I=1 ,MAX(2,NQC)) 

DO 10 N=1, NGRID 

READ(IUNIT) ((((Q(J,K,L,NV),J=1 ,JDIM(N)),K=1 ,KDIM(N)),L=1 ,LDIM(N)),NV=1 ,NQ) 
10 CONTINUE 


Special Notes 

(1) The number of grids is given by NGRID. The grid dimensions are given by 
JDIM, KDIM, and LDIM for the J, K, and L directions, respectively. 

(2) Single grid format for all of the above files is recovered by omitting the first 
READ statement and using NGRID=1 for the rest of the READ statments. 

(3) XYZ files without IBLANKS are also acceptable and are read by omitting the 
READ on the IBLANK array. 

(4) The solution file can be in PLOT3D or OVERFLOW format. NQ is the number 
of Q variables. The first five Q variables contain the density, X, Y, and Z 
momentum and total energy per unit volume, respectively. The sixth Q variable 
in an OVERFLOW solution file contains the ratio of specific heats. All other Q 
variables, if any, are disregarded. If a PLOT3D solution file is used, the ratio 
of specific heats is assumed to be constant everywhere and equal to GAMINF 
in the input parameters file. 

(5) The NQC,IGAM,HTINF,HT1,HT2, and RGAS parameters are disregarded. 

(6) Solution files for earlier versions of OVERFLOW are also acceptable. These 
may not contain the BETA and TINF parameters. 
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quantities integration ( OVERINT module). Instructions on how to use this software package are described in 
this user’s manual. Equations used in the flow coefficients calculation are given in Appendix A. 
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